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ABSTRACT 

Ilk this study ue develop flight director logic for flight path and 
airspeed control of a pouered-llft STOU aircraft in the approachr 
transition* and landing configurations. Tuo methods for flight director 
design are investigated. The first method is based on the Optimal 
Control node! (OCM) of the pilot. The second method* proposed here* uses 
a fixed dynamic model of the pilot In a state space formulation similar 
to that of the OCM* and Includes a pilot uork-load metric. 

Several design examples are presented with various aircraft* sensor* 
and control configurations. These examples shou the strong impact of 
throttle effectiveness on the performance and pilot uork-load associated 
ulth manual control of powered-lift aircraft during approach. Improved 
performance and reduced pilot work-load can be achieved by using dlrect- 
lift-control to increase throttle ef fectlveiiess. 

Flight path regulation In the presence of turbulence* glide-slope 
capture* and acceptable response to horizontal uindshear* can be 
achieved* equally well* by both methods. The second design method 
provides Improved control logic for manual control tasks by reducing the ? 

pilot's uork-load. It Is also easier to use and validate, and it can 
accommodate any linear model of the pilot. | 
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Chapter 1 
INTROOUCTIOH 


1.1 BACKGROUND MID RUR^aSI 

Manual control of Short Take-Off and Landing (STOL) aircraft Is 
generally more difficult than control of current Conventional Take-Off 
and Landing (CTOL) aircraft. Some of the difficulties are Inherent t 
landing on a short field with acceptable sink rates requires more 
preci:;e flight path and airspeed control* uhlie stability and control 
effectiveness are reduced at lou speeds. 

The high lift coefficients necessary for STOL operation are generated 
by increased uing circulation and by thrust vectoring that results from 
blouing engine exhaust over specially designed tral 1 ing-edge flaps, 
flight tests have shown t1*16] that pouered-Iift aircraft are 
characterised by : (1) sluggish flight path response to attitude 
changes* (2) operation on the backside of the drag curve (i.e.» uhere 
drag increases when airspeed decreases)* and (3) large changes In lift 
and drag with engine pouer setting. Doth attitude and power changes 
produce significant changes in flight path angle and airspeed* which 
makes manual control difficult. In the approach and landing 
configurations* pitch attitude is used primarily for airspeed control* 
and thrust modulation is used primarily for flight path control. 
Exactly the opposite technique is used in the cruise configuration where 
the effect of controls is restored to that of conventional aircraft. 
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Most pouered-lift aircraft incorporate a pitch-rate command/pitch- 
attitude hold augmentation system that effectively eliminates attitude 
control deficiencies [2»16]. The remaining problems* namely handling 
quality deficiencies associated with flight path and airspeed control* 
can be solved either by additional augmentation or by using a flight 
director. 

The purpose of this study is to develop flight director logic for 
flight path and airspeed control of a pouered-lift STOL aircraft in the 
the approach* transition* and landing cor figurations. 

The function of a flight director is to process aircraft sensor data 
and to display to the pilot appropriate pitch attitude Cor stabilator) 
and throttle commands. T‘ e flight director performs the cross-coupling 
and equalization that the pilot uould otheruise have to provide. A 
flight director is less expensive than an autopilot* and provides the 
pilot uith the training necessary for manual flight uith conventional 
instruments. For this latter purpose* the usual status information 
(i.e.* airspeed* radar altitude* glide-slope deviation etc.) are also 
displayed so that flight director failures can be detected* and the 
landing continued or aborted safely. 

For the flight director to be effective it must provide satisfactory 
performance and produce acceptable pilot work-load. First* flight 
director control commands* if followed exactly by an ideal pilot* should 
result in the desired aircraft response. Second* the displayed commands 
should be coroatible uith human operator capabilities so that the pilot 
can track and execute the commands with acceptable levels of mental and 
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physical efforts. The last requirement Implies knouledge about the human 
operator characteristics, uhlch Is the subject of manual control theory. 
Our knouledge about the pilot as a dynamic control component Is far from 
being complete or even extensive, and notions such as ^pilot uork-load^ 
still resist rigorous definition and treatment. Nevertheless, 
mathematical models of the pilot have proven useful in improving 
pilot/vehicle integration. 

Tuo sucii models are relevant to this study. The crossover model t4j 
Is the best knoun pilot-model for single-input single-output tracking 
tasks. Host of the research about pilot preference in control led-element 
dynamics has been done ulth this model [4, 5]. Several (single -input) 
flight directors were designed according to the results of this research 
which clearly Indicated pilot preference for plant dynamics having k/s 
like characteristics. The second model, known as the Optimal Control 
Model (OCM) of the pilot 161, is inherently capable of treating multi- 
variable systems, and Is used extensively In this study. 

The first flight director design method investigated In this study is 
based directly on the OCM. This method has been used before !9l, but the 
control configuration investigated there, a longitudinal hover task, is 
relatively uncoupled, and thus places less demand on the pilot than the 
STOL landing task investigated here. A second design method is proposed 
here, which Includes consideration of the pilot work-load, and thus it 
should produce Improved flight directors for manual control tasks. 
Several design examples are given to show the relative merits of the tuo 
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Manual control thaory is ravicuad in Chapter II. The concept* 

N 

structure* and special features of the OCd and crossover models are 
described. The aircraft and the control task are described in Chapter 
III. The Quiet Short-Haul Research Aircraft CQSRA) and its flight 
control system are presented* followed by a definition of thr landing 
task and the associated aircraft configurations. The flight director 
design requirements* including design constraints and performance 
criteria* are given in the last section. 


Two flight director design methods are presented in Chapter IV. The 
design procedure for each method is explained and the two methods are 
compared. 


The pilot/vehicle mathematical model is developed in Chapter V. In 
Chapter VI, both design methods are used in design examples with various 
aircraft* sensor* and control configurations. Conclusions and 
recommendations for future research are given in the last chapter. Model 
parameters for the QSRA are included in Appendix A. 
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Chapter II 

REVIEW or MANUAL CONTROL THEORY 

Until 1965 most of the research in the field of manual control uas 
devoted to understanding the characteristics of the human operator as 
the controller of a single-input single-output CSISO) tracking task. The 
result was a set of quasi-1 inear models that predict human behavior 
quite uell for these simple but important tasks. An excellent summary of 
this uork can be found in the report by McRuer et al . M). Tuo 
approaches have been used to extend manual control theory to the multi- 
variable case. The first, using classical multiloop control theory, 
relies heavily on judgments concerning the closed-loop system structure, 
and consequently is difficult to use in a systematic fashion. The 
second, knoun as the Optimal Control Model of the pilot (61. uses state 
space methods and optimal control theory. and is inherently capable of 
treating mu 1 1 i -var i abl e cases. 

2.1 .a.ASSlCfVL MANUAL CONTROL 

Classical manual control theory has been applied most successfully to 
time-invariant SISO compensatory tracking tasks. The model used to 
describe the pilot's behavior in these tasks is knoun as the crossover 
m odel |4J. The two essential elements of this model are (I) a 
describing-functi''.) analytical form, and (2) a set of adjustment rules 
which specify how to set the parameters of the model. 
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gure 1: Block Diagram of a SISO Tracking Task 
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A typical SISO laboratory tracking task is Illustrated in Figure U 
The pilot acts to minimize the displayed errors e> betueen the desired 
inputs i» and the aircraft response y. Tht actuator* airframe* sensor and 
display dynamics are all included in the control led-element dynamics 
represented by the transfer function YcCs). The pilot is represented by 
the quasi-1 inear describing function The remnant represents the 

non-linear part of the pilot behavior. The input forcing function* i* is 
modeled as a stationary random signal uith Gaussian distribution* rms 
level (7{* and banduith 

The crossover model relates the form of the pilot equalization* Yp* 
to the control 1 ed-el ement dynamics* Yc* by the equation : 

-j(lT 

Uce 

Yioltjwl = Yp(jw)Yc(jw) = near «c (2.1) 

3U 


where 

Y*ol is the system open-loop transfer function 

Uc is the crossover frequency (i.e.»uhere lYsoll^O dB.) 

T is the effective time delay representing transport lag 
and high frequency neuromuscular dynamics. 

Eqn.(2.1) has as its basis a large body of experimental evidence 
[^*51. Data from these experiments have shown that for a wide range of 
control 1 ed-el ement dynamics such as Ye(s) = Ke* Kq/s* and Kp/s^* the 
above relationship can be sat sfied uith a pilot describing function of 
the form : 

(Tijw+1) -j«T 

Yp(jw) = Kp e (2.2) 

(Tijw+l) 
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where 

Kp is the pilot gain 

N 

Ti is the lead time constant 
T; is the lag time constant 
T is the effective time delay 

The crossover model predicts that* given the control 1 ed-el ement 

dynamics Ye* the pilot will adapt his behavior in the following way : 

1. The lead and lag equalization* Ti and T;* are adjusted by the 

pilot to achieve a -20 dB/decade slope in the system open-loop 

amplitude response ( lYsoil ) near or below the crossover 

frequency «c- 

2. The pilot gain* Kp* is adjusted to locate the crossover frequency 
Ue=KpKc above the input forcing function bandwidth Wi so as to 
minimize tracking errors. 

3. The pilot effective time delay* r* is adjusted to provide 

adequate phase margin 

For the crossover model the phase margin is given by : 
i»iB = r/2 - T«c 

Experiments have shown that the operator tends to minimize the mean- 
square tracking error M*51. For a rectangular input spectrum of 
bandwidth and rms level ci* and under favorable tracking conditions 
(«c > Wi)* the rms tracking error is given by : 
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To decrease the rms tracking error the pilot tends to increase the 
crossover frequency Mq by increasing his gain Kp. However the pilot 
time del ay » reduces the phase margin and usually limits Cfe to less 
then 10 rad/sec. To increase the phase margin the pilot can decrease his 
effective time delay by concentrating more on the task and increasing 
his neuromuscular tension. This represents an increased work load* and 
in any case the pilot's effective time delay cannot be reduced below 
some physiological limit, usually about 0.1 sec. Uhen the pilot has to 
generate a low frequency lead, a larger time delay is incurred. 

A substantial amount of research has been devoted to finding the 
relationship between the parameters of the crossover model and 
subjective pilot ratings used to describe flying qualities. Good 
correlation between these parameters would permit use of the crossover 
model to design for desirable aircraft handling characteristics. The 
results of this research have indicated that the best pilot ratings are 
obtained when the pilot equalization is minimal (i.e., no lead or lag is 
required). Only moderate degradations in ratings appear when lag or 
small values of lead are required. If larger values of lead equalization 
are required the pilot ratings degenerate rapidly [41. 

2.2 oB um . ngoLi. ql mi 

The Optimal Control Model (OCM) of the pilot is based on the 
assumption that a well-trained pilot has the ability to estimate the 
state of the aircraft, and to produce control actions that minimize the 
integral -square output errors with a constraint on integral-square 

The main differences between the OCil and other 
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models of the human operator are the extensive use of state space 
concepts* the methods used to represent human limitations* and the use 
of modern control theory to compensate optimally for these limitations. 
The essential features of this model and the method of application are 
reviewed briefly in this section. A detailed description of the OCM 
concept can be found in the paper by Kleinman* Baron and Levinson (61. 

2.2.1 Model Description 

A birck diagram of the closed-loop pilot-vehicle system is shown in 
Figure 2 [9]. Since the OCtI method can handle multi-input multi-output 
CniflO) systems* all the variables shown in the figure are vector 
quantities. 

The displayed variables are assumed to be corrupted by an 
''observation noise" that accounts for the pilot's limitations in 
perceptual resolution and attention-sharing capacity. A Kalman filter is 
used to represent the pilot's ability to estimate the current state of 
the system. An optimal regulator is then synthesized to produce a set 
of control commands that minimize a quadratic performance index. The 
time delay element accounts for the pilot's central processing 
limitation. Just as observation noise is used to account for imperfect 
human perception* a "motor noise" is introduced to corrupt the optimal 
commands. This motor noise accounts for the pilot's inability to 
generate noise-free control actions and his imperfect knowledge of his 
own outputs. The last element of the OCM represents the bandwidth 
constraint of the neuromuscular system. 
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Figure 2: Block Diagram of the OCfl 
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Human eharacteristics# as incorporated in the OCn, oan be divided 
into four categories : 

\ 

1. Tr.sk formulation. 

2. Perceptual characteristics. 

3. Central processing characteristics. 

4. Neuromuscular actuation characteristics. 


2.2.2 Task formulation 

The pilot is assumed to adopt a strategy that minimizes e quadratio 
performance index of the form •- 


J 


» 



(y*Ayy + u^u)dt 


(2.3) 


uhere Ay and B are weighting matrices associated with the displayed 
outputs y and the controls u. The elements of Ay relate to the 
performance objectives of the particular control task (e.g.» glide-slope 
regulation)* while those of B reflect the control utilization 
constraints that the pilot or the aircraft systems may impose. This 
strategy is based on an extension of SISO laboratory tracking tasks that 
showed that the operator tends to minimize rms output errors while using 
specified control energy [4*6]. The OCM formulation* using quadratio 
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performance indices, has been validated for both StSO tasks and 

Hulti-Input Hul ti-Output (MIMO) tasks 

2.2.3 Perceptual Characteristics. 

The perceptual characteristics are modeled by t 
y = Hx 

yp = y + »)o 

where x is the state vector, y is the display vector, yp is the vector 
of perceived variables, Dq is the observation noise vector, end H is the 
output distribution matrix. In formulating the display vector, it is 
assumed that the operator can extract both position and rate information 
from a single moving cockpit indicator [6]. 

The statistical properties of the observation noise are determined 

by: 

1. Single observation. 

Studies of controller remnant have shown that the variance Qn of 
each white observation noise Doi is proportional to the variance 
of the associated perceived variable y{. A normalized noise- 
to-signal (N/S) ratio of -20 dB (i.e., correlation time of 0.01 
sec) is typically found in laboratory tracking tasks t7l. 

2. Task interference. 

In more complex situations, where several indicators are scanned 
by the pilot, the noise level is increased to account for the 
higher work-load. The scanning process appears to be highly 
complex and for most design applications a simplified form of 
task interference is used 111) : 
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fH = pQ/f\ 

Mhtre po is the besio noisc-to-signal ratio used for a single 

\ 

Indicatorf f{ is the fraction of attention devoted to indicator { 
(If{=l.)i and p{ Is the noise-to-signa) ratio associated uith 
indicator i. If attention is equally divided among all 

instruments* the noise to signal ratio uiU scale linearly uith 
the number of indicators. Finally* it is assumed that the 
perception of the rate of change of the displayed variables is 
realized without additional attentional demand* and hence* f{ 
does not have to be decreased on this account. 

2.2.4 Central Processing Characteristics 

Central processing characteristics include the pilot's equalization 
and a time delay. The equalization element is the most important part of 
the model as it represents the pilot's ability to adapt his behavior to 
the dynamics of the aircraft. The pilot's equalization is modeled by a 
Kalman filter followed by a linear quadratic regulator. This implies 
that tl<e pilot has an internal knowledge of both the dynamics af the 
system and the statistics of the disturbances. 

A single time delay is used to represent the accumulated transport 
time lag associated uith the central nervous system (visual* central 
processing* and neuromuscular transmission). Typical values are 
0.15>0.25 sec [5*6]. 
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2.2.5 MturomuscuUr Actuatlop Character istiaa 

The OCtI includes neuromuscular dynamics that account for 
physiological bandwidth limitations, and motor noise that accounts for 
the pilot's imperfect knowledge of his own outputs. Kotor noise is also 
used to represent the pilot's inability to precisely produce the desired 
motions. 


Neuromuscular dynamics are represented by a first order tag. 
Typically the neuromuscular time constant Tn is in the range 0 . 1 - 0.6 seo 
(41. The neuromuscular dynamics can be modeled explicitly as part of 
the control led-elcment dynamics. Alternatively, a control rate can be 
included in the performance index (eqn.(2.3)). The result is identical, 
but the explicit model eliminates the need to iterate the control rate 
weighting to achieve the desired value for Tn 

The motor noise is assumed to have the same form as the observation 
noise: the variance Rh of each white motor noise *)«{ is proportional to 
the variance of the associated control U{. tn the experiments 
reported in [81. a normalized noise-to-signal ratio of -25 dB Ci.e.» 
correlation time of 0.005 sec) was typically obtained by model matching 
anal ysi s. 


15 - 


OWGIIVAL' PAG£ Is 

OF POOR Ol’ALITY 


Chapter III 

AIRCRAFT AND TASK DESCRIPTION 


3. 1 AIRCRAFT DESCRIPTION 

The Quiet Short-Haul Research Aircraft (QSRA) is currently being used 
by the NASA Ames Research Center for terminal area» IoM-speed> pouered- 
lift flight research. The QSRA is a deHavilland C-8A Buffalo* modified 
by the Boeing Commercial Airplane Company to NASA specifications. It 
has a neu uing and nacelles and four AVCQ-Lycoming YF-102 engines 
mounted on top of the wing. The engines provide pouered-lift by 
deflecting exhaust gases over four specially contoured trai I ing-edge 
flaps referred to as the upper surface blowing (USD) flaps. Lift is the 
sum of the usual uing lift* the normal component of the thrust vector 
that results from flow turning, and the aerodynamic supercirculation 
created by the engine exhaust flow over the wing's upper surface. The 
aircraft has a maximum gross weight of 55i000 lb and a uing area of only 
6i)0 square feet. Maximum trimmed lift coefficients in excess of 10 have 
been demonstrated in flight. STOL landings are made at a lift 
coefficient of 5.5 which results in landing speeds of only 65-75 knots, 
while still maintaining a substantial stall margin {21. 
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3.2 FLIGHT CQNTRPL STSTEfl 

The airplane's control surfaces are Identified in Figure 3. The 
primary flight controls consist of a stabilator* rudder» ailerons* 
spoilers and engine power. Secondary controls consist of four USB 
trai I ing~edge flaps and two conventional outboard trai I ing-edge flaps. 
A three-axis stability augmentation system is provided. 

The use of the longitudinal controls is described below t 

1. The USB flaps are used for configuration control. The USB flap 

deflection is usually 0-10 degrees for takeoff* 0 degrees for 
cruise, 30 degrees for initial approach and go-around* and 50 
degrees for landing. USB flap deflection frnm 0 to 30 degrees 
is controlled by a handle in the overhead console. Flap 
deflections between 30 and 66 degrees are selected by a switch on 
the throttle which commands flap motion at a rate of 10 deg/sec. 
To enhance spanuise wing loading, the conventional outboard 

trai 1 ing-edge flaps are deflected to 59 degrees for all 

configurations except cruise. 

2. Pitch control is achieved by stabilator deflection. Uhen 
engaged, the pitch control augmentation system (PCAS) provides 
very tight rate command/attitude hold control. This type of 
pitch axis augmentation system is widely used on this class of 
aircraft and has the effect of increasing the short-period 
frequency and damping the phugoid mode. Attitude stabilization 
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Is rtquired In tht approsch configuration sine* changes In thrust 
or USB flap setting produce sizeable pitching moments. 

3. Thrust is controlled by a throttle in the usual usy. except that 
the orientation of the incremental thrust vector Is strongly 
influenced by the USB flap setting as well as by the nominal 
pouer setting. 

4. Oirect-1 if t-control (DLC) is achieved by use of symmetric spoiler 
deflections commanded by an electric interconnect from the 
throttle. In this mode the spoilers are first biased up to a 
setting of 13 degrees and move from that nominal position in 
response to throttle movements. To ensure that authority is 
maintained, a washout is norm.’tlly included. The higher bandwidth 
of the spoiler actuators and the increased effective control 
sensitivity result in faster flight-path response than can be 
obtained through the use of thrust modulation alone. 
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3.5 11^ Mimusm 

th« tasks invvstigstsd this study sr* ths initial approach* 
transition, and landing phasas of a poMtrad-*lift STOl aircraft. A 
typical approach and landing pattarn for th* ^SRA is shown bslow i 




C/ST0L CT0L 



DISTANCE 


I 


Thraa aircraft con f » gurat i on.'f 1 i gh t conditions ara used to represent 
tins task s 

I. CTi'.L - A str.-xi(3ht and level approach at an altitude of 1000 feat* 
and airspeed of 1.10 knots. The USB flaps are ratraeted. 

C^SJOl - In tins triinsilion configuration the USD flaps art 
deflev'ted .10 degrees. The aircrait is on a six degrees glida- 
s 1 opa . at 'TO K t . 

3. STOL - for the landing configuration tha USD flaps are daflaotad 
to SO degraes. The aircraft is on tha s.ima glida-slopa* at 70 kt. 
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The aircraft stability derivatives in the three configurations are 
given in Appendix A. 


3.4 DESIGN REQUIREMENTS 

The objective of this investigation is to develop flight director 
logic for flight path and airspeed control that uill provide good 
performance with iou pilot work-load in the approach and landing task 
described above. Special implementation considerationSf such as 
simplifying the control laws for real-time operation* and display 
design, are not addressed. 


3.4.1 Design Constraints 

To conform to the standard utilization of controls, the flight 
director should provide only throttle and pitch attitude commands. 
(Pitch attitude command is used as a control input rather than 
stabilator deflection; this simplifies the logic and is a good 
approximation because of the high bandwidth of the pitch control 
augmentation system described earlier). The USB flaps and the 
conventional outboard flaps are held fixed at the appropriate setting 
for each configuration. The spoilers can be used, if necessary, only in 
the DLC mode in which they are linked to the throttle and so can be 
modulated without increasing the pilot's work-load. These design 
constraints also allow the flight director to be incorporated without 
any modification to the existing control mechanization. 
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Quantitative performance criteria are stated in terms of flight-path 
and control deviations from trim in the presence of random wind 
disturbances* and in terms of the control authority required to capture 
the glide-slope. These criteria are based on flight test data obtained 
with an aircraft similar to the QSRA 1161. 

1. Disturbance rejection. 

Horizontal and vertical uind gusts are approximated by shaping 
filters driven by white noise with rms outputs of 2.3 ft/sec. 
The average output and control deviations (in the rms sense) 
associated with these disturbances were chosen not to exceed the 
values shown in Table 1. 

To ensure satisfactory transient responses to uind gusts* 
additional criteria are used. The maximum values of flight path 
and control deviations following two-sigma horizontal and 
vertical wind gust impulses should not exceed those shown in 
Table 1. 

2. Gi ide-slope capture. 

The capture of the glide-slope is represented by the aircraft 
response to an initial offset of 4 knots in airspeed or 20 feet 
in vertical position. Control deflections under these conditions 
were chosen not to exceed the maximum values specified in Table 
1> and residual flight path deviations after 20 seconds should be 
less then lOJi of the maximum values shown in the table. (The 
aircraft response to an initial offset in airspeed can also be 
regarded as the aircraft response to horizontal uindsheer). 
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Use of the flight director should reduce the pilot's uork-load. A 
quantitative criterion is not specified. A metric for assessing pilot 
uork-load is presented in Chapter IV. 


TABLE 1 

Flight-Path Regulation and Control Authority Criteria 


variable 

rms 

value 

maximum 

value 

u. 

[ft/sec J 

2.25 

6.75 

d 

[ftl 

6.7 

20. 

• 

d 

1 ft/sec 1 

2.2 

6.5 

8t 

l%] 

3 

6 

8 

Ideg 1 

2 

4 
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CHapttr IV 

FLIGHT UIRECTOR DESIGN METHODS 

Tho methods for flight director control tau design for • STOL 
aircraft are presented in this chapter. The first methods proposed by 
Levinson [10]» uses the OCM concept to predict the equalization 
characteristics of a uell-trained pilot performing the control task 
using only conventional cockpit instruments. These equalization 
characteristics, represented in the OCM formulation by a Kalman- 
fi1ter/LQ6 regulator, are then used as the control laws for the flight 
director . 

The second method, proposed here, uses a vehicle model which includes 
f 1 jqht director represented by a Kalman~f i 1 ter/LQ regulator, and a 
fixed dynamic model of the pilot. This method uas developed to overcome 
reservations about the first design method, in particular the lack of 
explicit design guidelines for reducing pilot work-load. 

The tuo design methods are described below and then compared to each 
other . 
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4.1 fllGKT DIRECTOR DISPLAY 

A typical flight director display is shown in figure 4. 

N 

The ''director" part of the display includes : 

1. Throttle command bar (a) 

2. Pitch attitude command bar (b) 

The "status" part of the display includes at least the fotlouing t 

1. Aircraft symbol Cc) 

2. Artificial horizon with pitch attitude scale (d) 

3. ILS box (g1 i de-*s1 ope and localizer) (e) 

4. Airspeed (f) 

5. Radar altitude (g) 

6. Sink-rate (h) 

The status information enables the pilot to detect flight director 
failures* and to continue the flight safely. 
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Figure 4: Flight Director Display 
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4.. DESIGN neiHgP 1 ^ 

Block-diagrams of the pilot/vehicle system for the design and 
implementation steps of method 1 are shoun in Figure 5. For the design 
step* the pilot is assumed to be using conventional cockpit Instruments 
(e.g.* airspeed* sink-rate* and glide slope deviation indicators). Using 
the OCn concept as outlined in Section 3.2* the pilot's equalization 
characteristics are represented by the Kalman-f i1 ter/regul ator block. 
The Kalman-f i 1 ter and regulator are then regarded as a compensator* and 
are implemented* either directly or preferably in a simplified form* as 
the flight director control laws. Thus configured* the flight director 
will provide the same throttle and pitch attitude commands as would a 
uel l-trained pilot, 

A point to note is that the pilot's time delay and neuromuscular 
dynamics are not included in the flight director control laws. These 
characteristics are assumed to be tasK'invariant and will be "supplied" 
by the pilot regardless of the form of the display compensation. On the 
other hand* including the time delay and neuromuscular dynamics in the 
model used for the design step* ensures that the resulting control laws 
will compensate "optimally" for these limitations. 
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AIRCRAFT 

X 

COCKPIT 

V 


INSTRUMENTS 




STEP 1 : DESIGN 



PILOT 


DISPLAY 


STEP 2 : implementation 


Figure Pitot/Vehicle Block-Diagrams for Design method I 
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4.2.1 Dtaiqn Proctdur# 

To satisfy both tht OCtI specifications and the flisl^t director design 
criteria the follouing procedure Is used : 

1. A state space model of the pllot/vehlolc system Is set up 
according to the OCn formulation (o.f. Chapter V for modeling 
details) . 

2. Appropriate values for the output and control weighting matrieesr 

Ay and of the cost function J (eqn.(2.3)) are selected to 

reflect the objectives of control task. This subject ulU be 
discu.;^ed belou. 

3. Initial values for the spectral density matrices^ Q and R (the 
pilot motor and observation noises) are determined. Given the 
wind gust statistics used in the model. ie>;pected values for the 
rms outputs and controls are assumedf based on aircraft 
characteristics and previous flight experience (e.g.* the flight 
test data reported in 1161). These expected values are used to 
calculate the elements of Q and R that satisfy the 0CI1 
requirements for p*1ot motor and observation noise-to-slgnal 
ratios (c.f. Section 2.2.3). 

4. The OCn is generated by computing the Kalman-f i 1 ter and LQ 
regulator gains. 
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5. Tht spectra) density matrices* Q and R* are adjusted* and step 4 
is repeated until the specified noise-to-aignal ratios are 
obtained. 


6. The output and centre) weighting matrices* Ay and B* are 
adjusted* and steps 4 and 5 are repeated as necessary unti) the 
f)ight director design criteria* set forth in Section 3.4* are 
met . 


Setting up the state space pi )ot/vehic)c mode) and adjusting the 
matrices Q and R to achieve the specified noise-to*signal ratios is 
re)ative)y straightforward. The essentia) probiem is therefore the 
proper se)ection of the weighting matrices of the the cost function. 

4.2.2 Select ion ai H e i ah LLn.a natri.gg-S 

In the design examples reported in [9] and ()0]» the weighting 
matrices were based on Bryson's ru)e (12). Thus the diagonal elements of 
Ay and B are equal to the inverse of the square of the maximum allowable 
output deviations and control usage respectively. In both of these 
design examples the weighting matrices were held fixed* i.e.* the last 
step of the design procedure presented above was not used. The design 
reported in [9] was evaluated in a piloted simulation which showed 
reasonable agreement between the rms outputs and controls predicted by 
the model and those obtained in the simulation. This evaluation also 
snowed a definite reduction in rms outputs (and a smaller reduction in 
rms controls) for the flight director configuration compared to the 
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conventional display configuration and an improvement, in pilot opinion 
ratings. For this study* using rms outputs and controls as the sola 
criterion was deemed insufficient* especially since a piloted simulation 
was not available to validate the design. Thus criteria based on 
c1osed“loop pilot/vehicle transient response were added (c.f. Section 
3.4). The procedure for adjusting the weighting matrices to meet the 
criteria is presented in Chapter VI. 


4.3 DESIGN HETH0D 11 

A block-diagram of the pilot/vehicle system for design method II is 
shown in Figure 6. The vehicle model includes a flight director 
represented by the Kalman-f i I ter/regul ator block. The pilot is 
represented by a fixed dynamic model* with unity gain equalization* that 
is assumed to represent the lower work-load and improved pilot opinion 
associated with the use of a flight director. This pilot mode! is an 
extension of SISO manual control experimental results. These experiments 
have shown (c.f. Section 2.1) that the best pilot opinion ratings are 
obtained when the pilot equalization characteristics are minimal (i.e.» 
unity gain equalization). If the flight director is well designed* the 
pilot's task is simply to transfer the two displayed commands to the two 
aircraft controls on a one-to-one basis* subject to his inherent time- 
delay and neuromuscular lag. Under these conditions it is assumed that 
the experimental results* although derived for a SISO task* are still 
useful design guidelines. 
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Figure 6: Pilot/Vehicle Block-Diagram for Design Method II 
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The pilot is thus represented by e unity gain equal ization^ folloued 
by a time*de1ay and a neuromuscul nr lag that are modeled in the same way 
as in method I. The OCM formulation is also used to model the pilot 
observation characteristics with the follouing changes : (1) the 
observation noise levels are louer since only two variables are 
displayed to the pilot and consequently less scanning is required (c.f. 
Section 2.2.3), and (2^ the 'ilot observation noise is treated as a 
process noise. Finally, sensor .^oise is introduced to corrupt the flight 
director inputs. (Sensor, or instrument noise was not included in the 
design model for method ! since it is negligible compared to the pilot 
observation noise). 

4.3.1 Work load Netric 

The pilot model used in method II is derived from a particular case 
(unity gain equalization) of the SISO crossover model (c.f. Section 
2.1). According to the crossover model the pilot equalization 
characteristics and the control 1 ed-el ement dynamics are related by the 
system open loop transfer function in eqn.(2.1). These experimental 
results indicate that the system open loop transfer function has a 
strong effect on the work-load associated with the control task. It is 
suggested here that the system open loop transfer functions will have a 
similar effect in tllMO systems. Experimental measurements of frequency 
response in Mino systems are quite difficult and consequently the form 
of such a relationship has not been determined yet. (The validation of 
this design method in a piloted simulation may provide the answer). 
Therefore a modified form of eqn.(2.1) is used as a work-load metric i 
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for tMch control ohanntl. th« gain from the control input to tho 
aircraft to the commanded control output of the pilot* ehould behave 
like k/t in the frequency range of intereat. Cross transfer funotions 
and phase characteristics are ignored. 

rxerpt for the additional uork-lond criterion* the design procedure 
for method II is the seme as that used in method 1 and is not repeated 
here . 

4.4 £QtJP.\!US0H ai IM BiSliia 

The second design method uas developed to overcome reservations and 
possible deficiencies in using the first method. These reservations 
arose from the fact that design method I i.s unusual since the situation 
being modeled in the design stage is not the one in which the design 
Mill be used. Specifically* the flight director control laws are based 
on a model of a pilot perfeming a high work-load task using conventional 
cockpit displays. The Kalman-f i 1 ter gains (and thus the flight director 
control laws) depend among other things on the pilot observation noise. 
The pilot observation noise levels postulated by the 0C11 are much higher 
than the sensor noises that are corrupting the flight director inputs. 
Thus, the Kal man-f i I ter is designed for a measurement noise level much 
higher than the one it will actually experience. Another deficiency is 
that the first method does not include any mechanism to ensure a lower 
pilot work-load other than that implicit in the concept itself (i.e.* a 
flight director control law that mimics the behavior of a well trained 


pilot IS used! . 
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A piloted simulation uould normally be required to evaluate a flight 
director before actual implementation. Houcver» this type of validation 
is not well suited to the several design cycles typically required In 
the development of advanced flight director systems. Instead* some 
method to permit a preliminary analytical validation is highly 
desirable. Since design method I is not based on the actual situation 
in uhich the design uill be used* a separate analytical validation is 
especially desirable. Such a validation procedure obviously requires on 
analytical model of the pilot* and the one used in design method II is a 
possible candidate. Design method 11 is thus more efficient as the same 
mathematical model can be used for both the design and the validation 
phases and unnecessary design iterations are eliminated. 

Another advantage of the second design method is that any (linear) 
model of the pilot can be used, whereas the first method is limited to 
the OCfl. If further research in manual control theory results in a 
better understanding of the human operator, the improved model can be 
incorporated in design method 11 with only feu modifications. 
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Chapter V 

AIRCRAFT-PILOT MODEL DEVELOPMEMT 

In this chapter the pilot/vehicle state space model is developed. Most 
of the chapter is devoted to the development of the pilot/vehicle model 
for design method I. The modifications required for design method 11 are 
presented in the last section. 

The pilot/vehicle model for design method I (c.f. Figure 2) includes 
the foUouing : 

1. A tuo degree-of-f reedom translation model of the airframe. 

2. A vertical and horizontal wind gust model. 

3. A first order thrust model. 

4. A fourth order model of the combined pilot time delap and 
neuromuscular systems. 


These elements represent the open-loop plant. The 
process is represented by a Kalman-f i 1 ter/regul ator 
formulation. 


pilot equalization 
in the usual LQC 
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5.1 AIRFRAME STATE EQUATIONS 

The linearized equations for longitudinal motion are t 

u “ XuU “ (gcos6|}}6 + + X^8s + 

- (XuCOs6o+XHsin8o)UM ~ (XuSindg-^XMCosSoIUM 


(5.1) 


u - ZuU ZmU - (gsin8o)9 (Uo-**Zq)q Zq8« Z^S. Z^pGap 

(5.2) 

- (ZuCOs9o-t-ZMSin6o)UM - (ZuSin9o-*’ZHCos9o)UM 


d = “U + Uo9 (5.3) 

8 = q (5.4) 


where 

u>u are horizontal and vertical velocity perturbations 
Um'Um are horizontal and vertical wind gusts 
q>8 are pitch rate and angle 

d is vertical perturbation from the glide-slope. 

Ga>Gap are stabilator and spoiler deflections 
Sa is engine speed 

X()>Z() are axial and vertical force derivatives 
Uo is x-axis initial velocity component 
8g is initial pitch angle 

Stability axes are used. The quantities u,Uh and X are defined positive 
forward. w.Wh and Z are positive downward; d is positive upward . 
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The pitching moment equation is not included in the above set since 
ue have used the fact that the QSRA aircraft is equipped uith a very 
tight pitch control augmentation system. This simplifies the aircraft 
control problem from three to two degrees of freedom. Consequently 
commanded pitch angle 8a is regarded as the control input instead of 

which can nou be eliminated on the implicit assumption 0 > 8«. To 

* 

eliminate the term Uo6 from the equations* d = -(u-Uq 8) uas used as a 
state variable instead of u. 

The following usual assumptions were made to further simplify the 
equations : 

• 6o small* thus cos8o-l and sinSo^^do* 

• Zq small compared to Uq. 

• X„sin0oUH etc. are negligible. 

• The stabilator contributions to X and Z can be neglected. 

Uith the above assumptions the airframe state equations are reduced 
to : 


u " X ijU “ X vid + X (S 0 X^pSgp (XMUo~g)8a “ X^jUm “ X^Wn C 5 . 5 ) 

• « • 

d ~ Z yU ~Z| 4 d ^ Z (S • ^ ^sp^sp ^ C Z 14 U 0 ~ g 8 0 ^ 8 n ~ 2!yUi4 "* Z^w^ C5.6) 


Note : In the following, the subscript a of 8« will be dropped. 
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S.2 TURBULENCE SHAPING FILTERS 

Uind gusts are modeled as first order Harkov* processes. The 
relationship between the gust parameters and aircraft geometry* airspeed 
and altitude was taken from Holley and Bryson M4I. 


5.2.1 Horigontal wind oust 

Um(s)/t)u(s) = 1/(TuS+1) 
or 


Uh = 



T»U 


Tu Tu 


where 

7)u = white noise with zero mean and 
spectral density Qu 
Qu = 2o*L/CiV 
Tu = L/C,V 

with 

0 = rms(UN) =2.3 ft/sec 
V = airspeed 

L = turbulence integral scale = Lah/(ho+h) 

La = 2000 ft 

ho = 2500 ft 

b = wing span 

B = b/2L 

( 1 + 30 / 2 )*^’ 

C, = 

1+3B 


(5.7) 
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5.2.2 Y.ffn.Lis.al H.ind sml (Hii) 

The simplified form for small B uss used to obtain a first order gust 

N 

mode' t 

Wh = + — (5.9) 

Tn T„ 

Hhere 

Dm = white noise with zero mean and 
spectral density Qw 
Qh = »*L/V 
Tm = L/2V 

5 . 3 thrust ru ?G£l 

In the flight regime of interest the thrust ic assumed to follow 
engine rotational speed linearly. The engine speed response to throttle 
position is modeled as a first order leg : 

• • St 

S, = + — (5.10) 

T. T. 

where 

St = throttle position, in units of 
equivalent engine rpm 
T» = engine time constant 
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5.4 DIBEgT nil CONTROL P1ECHANI2ATIQN 

In the OLC mode the spoilers are linked electricaly to the throttl 
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5.5 E.1LQ.T UnmUAT AM HEgKQnVSCVLAS 0TI<ftmCS STSTLn 

Thest two tlements are connected as follou : 



noi$d 


uhere 

8* = desired control (regulator output) 
Sd = delayed control 
8 = control input (to the aircraft) 
t>bi - pilot motor noise 


The time-delay element is modeled by a Pade approximation i 
6d(s) -(s-2/t) 

= ( 5 . 11 ) 

8*(s) (s+2/t) 

The neuromuscular system is approximated by a first orJer lag t 
8(s) 1/Tn 

= ( 5 . 12 ) 

8d<S)+’)mCs) S+l/Tn 
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The transfer function betueen the desired control end the actual 
control input is obtained by combining the tuo elements : 


-I/T„(s-2/t) 1/Tn 

S(s) = 8*(s) T U„(s) 

Cs+l/TnHs+2/T) (s-M/Tn) 


-l/Tn(s-2/T)6*Cs> ♦ 1/Tn(s+2/T)fl«(s) 

(5.14) 

S* + (2/T-H/Tn)« 2/(rTn) 


There are tuo such systems* one for the throttle channel (8«) and one 
for the pitch attitude control channel (6). Transforming eqn. (5.14) 
into a state space form with the appropriate parameters for the throttle 
channel ue get : 


« 

St 'I r -C2/Tt+1/Tnt) 
•StJ -2/(TtT„t) 
r -l-'Tnt 

* 

L 2/(TtT„t) 



St* r 1/T„t 

L 2/(rtTnt) 



And correspondingly for the pitch control channel i 



r -(2/Tt+1/Tni) 
L -2/(r*Tn#) 



♦ 


r "l/Tni 


r l/Tne 

l^lae 

2/(t,T„,) 

] * 

■ 2/(rfTne) 

] 


iS.15) 


(5. 16) 
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8tb = throttle channel internal state 

\ 

9b = pitch channel internal state 
Ti = throttle channel time-delay 
r« = pitch channel time-delay 

Tnt = throttle channel neuromuscular time constant 
Tnt = pitch channel neuromuscular time constant 

5.6 PILOT/VEHICU 8.YSTEM FOR DLSIMI flLrHDQ I 

Combining the various elements, the following dynamic system is 
obtained (c.f., top of Figure 5) : 

X = Fx + Gu + Ftj 
y = Hx + Uo 

where 

Ihe state vector x is : 


u 

(ft/sec I 

Horizontal inertial velocity 

• 

d 

(ft/sec 1 

Sink rate 

d 

Ift) 

Vertical position 

s. 

IJCl 

Engine speed 

Si 

ly,] 

Throttle position 

S tb 

15',/sec ) 

Internal throttle state 

9 

Idegl 

Pitch attitude 

Ob 

( deg/sec ) 

Internal pitch attitude state 

Um 

1 f t/sec 1 

Horizontal wind gust 

Um 

t ft/sec ] 

Vertical wind gust 
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The control vector u is t 

(Xl Throttle commend (regulator output) 

6* tdeg] Pitch command (regulator output) 

The process noise vector t) is : 

Turbulence 
Pilot motor noise 

The output/observation vector y is chosen to be : 


U» 

1 f t/sec 1 

Airspeed 

d 

If tl 

Vertical position 

• 

d 

1 f t/sec ) 

Sink rate 

d 

1 f t/sec‘ 1 

Vertical acceleration 

u 

1 f t/sec* 1 

horizontal acceleratii 


The associated observation noise vector is t 

* ** • 

TJou' ’lod' ^od’ ^od» ^ou- 


The actual number of outputs used in a given design may vary. 
See Chapter VI for details. 
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The state matrix is : 


Xy 

Xm 

0 



0 

XhUo-S 

0 

"Xy 

“Xm • 

-Zu 

-Zm 

0 

-Z. 

"kZ^p 

0 

-Z^Uo+OYo 

0 

Zu 

Zh 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1/T. 

1/T. 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2/Tt+l/Tnt 

1 

0 

0 

0 

0 

0 

0 

0 

0 

-2/(TtT„t) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-2/Ta+1/Tna 

1 

0 

0 

0 

0 

0 

0 

0 

0 

**2/(TtT n#) 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-1/Tu 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-UTh ■ 


The 

output 

distribution 

matrix 

is : 




1 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

Xu 

X w 

0 

X. 

kX jip 

0 

XmUo-9 

0 

-Xu 

-Xm 

-Zu 

"2(4 

0 

-Z. 

“kZsp 0 

-ZwUo+gyo 

0 

Zu 

Zh 


: k = 0 if DLC is not used. 
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The control and process noise distribution matrices are : 


0 

0 1 


P 0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


0 

Q 

0 

0 

0 

0 


0 

0 

0 

0 

“ 1 n t 

0 

r = . 

0 

0 

nt 

0 

2/(t iT nf) 

0 


0 

0 

2>^(TiT„i) 

0 

0 

-VTna 


0 

0 

0 

I^T na 

0 

2/(T*Tn9) 


0 

0 

0 

2/(rBTn#) 

0 

0 


1/Tu 

0 

0 

0 

0 

0 


0 

VT 

K 0 

0 


5.7 PllOT/VEHICLE SYSTEM FOR DESIGN METHOD JJ. 

The pilot/vehicle system for design method II is the same as the one 
for design method I except : 

1. The aircraft outputs, instead of being displayed to the pilot, 
are used as input to the flight director. 

2. Only two variables are displayed to the pilot : Si* and $* . The 

pilot observation noise is now considered as a process noise. 
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Referring to Figure 6* the follouing system is obtained : 
X = Fx + Gu + Ft} 
y = Hx + D, 

where 

The process noise vector t> is : 

turbulence 

pilot motor noise 

pilot observation noise 

The measurement (sensor) noise vector is : 

* ** 4 

T)su> ‘JJsd» I'lsd' flsd» flsu* 


Only the process noise distribution matrix is different : 


r 


r = 


0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 1/Tnt 

0 0 2/(TtTnt) 

0 0 0 

0 0 0 

1/Tu 0 0 

0 1/Tm 0 


0 

0 

0 

0 

0 

0 

na 

2/(T9Tn9 ) 
0 
0 


0 

0 

0 

0 

-VTnt 

2/(rtTnt) 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

-t/Tn* 

2/(T9Tn« 

0 

0 


) 
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Chapter VI 
DESIGN EXAMPLES 

Several flight director design examples are presented in this 
chapter. The first section describes the software tools used to 
generate these examples. The main design effort uas a trade-off study 

on the number of aircraft sensors and the amount of control authority 
required. The STOL configuration uas used in this study to gain insight 
about the tuo flight director design methods presented in Chapter IV. A 
practical procedure for satisfying the design criteria by iteration on 
the cost function weighting matrices was also establi ihed. 

The trade-off study examples are presented in the follouingi order : 

1. Three sensors configuration - Baseline Characteristics t all 
performance criteria except one are met ( insufficient throttle 
effectiveness is the limiting factor); pilot work-load level 
judged too high. 

2. Five sensors configuration - Significant reduction in pilot work- 
load; small improvement in performance. 

3. DLC configuration with three sensors - Significant improvement in 
performance due to increased throttle effectiveness; also pilot 
work-load is similar to that of the five sensors configuration 
but only three sensors were required. 
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The last section of this chapter includes a simitar design example for 
the C/STOL configuration. 

6.1 S0FT14ARE TOOLS 

A computer program, "OPTSYS'', developed by the Department of 
Aeronautics and Astronautics at Stanford University, uas used to 
determine the Kalman-'f i 1 ter and regulator gains 1131. The program uas 
modified to calculate the pilot observation and motor N/S ratios, and to 
predict the appropriate values of the spectral density matrices Q and R 
that would result in the N/S ratios specified by the OCM method. Uith 
this modification, the specified N/S ratios could usually be obtained 
uith only tuo iterations. 

The original program provided a printout of rms states, control s, and 
outputs uhich uere checked against the performance criteria (c.f. 
Section 3.4.2). The program uas further modified to provide ; 

1. An output disk file containing the closed-loop state matrix 
computed by OPTSYS and the control distribut'ion matrix. These 
matrices uere used as input to a time response program written by 
the author. The time response printouts and plots uere then 
compared to the performance criteria on wind gust rejection and 
glide-slope capture given in Section 3.4.2. 

2. A printout and output disk file containing the system open-loop 
( Y»ol 1 frequency response from 0.1 to 100 rad/sec. These data 
uere used to assess the pilot uork-load and satisfy the 
additional criterion proposed for method II in Section 4.2.1. 
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6.2 PILOT PARAMETERS AND SENSOR NOISE 

6.2.1 Pilot p . oran<^ t ers 

The pilot time delay and neuromjscul ar lags used are t 
Tt = 0.15 sec Te = 0.20 see 

Tni = 0.20 sec Tnt ~ 0.25 sec 

6.2.2 Sensor Noise 

The values chosen for the elements of R (for design method II only) 
are based on conservative estimates of typical aircraft sensor noise 
levels. The noise sources are independents so R is diagonal and each 
element has the form ; 

Rii = 2T,Oi* 

uhere cr; is the rms noise level for sensor i* and the correlation time 
Ts is chosen to be smaller than the sample interval for typical real- 
time Kal man-f i 1 ter implementation. 

The values used are : 

Tm - 0.025 sec 
= 0.75 ft/sec 
Cd = 1. ft 
od = 0.5 ft/sec 
Od = 0.2 ft/sec* 

Ou = 0.2 ft/sec2 


51 


6.3 TRADE-OFF STUDY 


ORIQTNAL PAGE fS 
OF POOR QUALITY 


6.3.1 Three Sensors Configuration 

Three sensors uere used for the initial design example t airspeed 

e 

(ua>/ vertical position (d)> and sink-rate (d). In design method I 
these variables are displayed to the pilot by three cockpit instruments 
(c.f. Figure 5). The basic pilot observation N/S ratio (-20 dB) Is 
multiplied by three to account for the scanning process^ and results fn 
an effective normalized N/S ratio of -15 dO. In design method II these 
three measurements are processed by the flight director* and only tuo 
variables are displayed to the pilot : throttle command (8t>» and pitch 

attitude command (6). The pilot observation noise level is therefore 
lower by one third (-17 dB). 

An initial set of cost function weighting matrices was obtained by 
Bryson's rule (c.f. Section 4.1.2)* based on the maximum errors 
specified in Table Llimits.. The resulting flight directors* using 
either method* satisfied both the statistical and deterministic wind 
gust rejection criteria. The glide slope capture and pilot work-load 
criteria were not satisfied* and the pitch attitude control was felt to 
be under-utilized. 

The adjustment procedure for the weighting matrices that evolved 
during this study is : 

1. General - The pitch attitude weighting was used as reference and 
was not modified. 
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2. Performance criteria - The airspeed error (u«) and vtriioal 

position error (d) ueightings uere increased in order to decrease 

N 

the maximum and residual errors on and d. The throttle 

ueighting uas increased to more nearly equalize the control 
utilization of the throttle and pitch attitude. The u» and d 

ueighting uere increased accordingly to keep the same ratio of 

outputs to controls ueighting. The sink rate ueighting had a much 
smaller effect and uas not modified. 

3. Mork-load criterion (method II only) - The effect of the 

ueighting matrices on the uork-load metric uas more complicated. 
No analytical procedure is known to a1 lou systematic shaping of 
the frequency response resulting from LQ6 design. Consequently 
trial and error adjustements of the ueightings uere made and the 
effects assessed. Taken individually* any ueighting except the 
one on sink rate affected mostly the lou frequency range of the 
system open-loop frequency response magnitude curve. 
Collectively* the previously described adjustments* required to 
improve the flight director performance* had a detrimental effect 
on the uork-load metric. That is* improving performance resulted 
in the frequency response magnitude curves becoming steeper and 
farther away from the required -20 dB/decade slope. The uork- 
load metric could be improved significantly only by increasing 
the sink rate ueighting or decreasing the sink rate measurement 
noise. These adjustments provided relatively good control of the 
frequency response magnitude curve slope for the pitch attitude 
channel, but insufficient control of the corresponding curve for 
the throttle channel. 
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The adjustment procedure was carried out separately for the tuo 
design methods. The end result was two sets of weighting matrices « one 
for method I> the other for method II. The difference between the tuo 
sets was not significant, and only one set (the weighting matrices as 
adjusted for method II) was used to generate the design examples 
presented below. 

Examples of flight director performance for the STOL configuration* 
using three sensors are presented in the f i rst column of Tables 2 to 6, 
and in Figures 7 to ID, 

As mentioned earlier, the wind gust rejection specifications were 
achieved. Table 2 shows the average outputs and controls in the presence 
of random wind gusts. Both methods produced similar results. Tables 3 
and A show the aircraft response to a two-sigma impulse in horizontal 
sr.d vertical wind gusts respectivly. A more detailed picture of the 
aircraft response is shown by the corresponding time history plots In 
Figures 7 and B. 

Glide slope capture following an initial offset in airspeed (4 kt) 
and vertical position (20 ft) are shown in Tables 5 and 6 respectively. 
The corresponding time history plots are shown in Figures 9 and 10. The 
glide slope capture criteria are met with one exception : following an 
Initial offset in airspeed, the vertical position residual error 
(2.18 ft) slightly exceeds the specified value (2.0 ft), while the 
maximum throttle excursion (6.30 %) is already above the 6. X limit. 
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As shown by the frequency response curves in Figure 11# the work'^loed 
criterion is not met either : whereas the pitch channel nagnitude curve 
slope (-18.5 dB/decade between 0.4 and 4.0 rad/seo) is quite acceptable* 
the coresponding slope for the throttle channel (-28.8 dB/decade) is 
much steeper then the desired -20 dB/decade. 

Although the performance criteria were not fully net* the main 
concern at this point was to improve the work-load metric. To this 
effect the use of additional sensors to measure higher derivatives of 
the aircraft outputs was investigated. 
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Mind Cust Rejection 


RMS Outputs end Controls 
STOL - nsthod 11 



outputs 

3 

5 

3 

design 


OLC ? 

no 

no 

yes 

criteria 

Ua 

1 ft/sec ] 

2.22 

2.20 

2.02 ! 

2.25 

d 

tftl 

5.00 

4.79 

3.07 

6.70 

• 

d 

1 ft/sec 1 

1.18 

1.13 

.974 

2.20 


\y.] 

2.29 

2.26 

1.59 

3.00 

e 

[degl 

.617 

.569 

.457 

2.00 



STOL 

• Method 1 




outputs 

3 

5 

3 

design 


DLC ? 

no 

no 

yes 

criteria 

u. 

I f t/-.ec 1 

2.24 

2.22 

2.04 

2.25 

d 

Iftl 

5.44 

5.08 

3.52 

6.70 

• 

d 

I ft/sec ] 

1.29 

1. 18 

1.13 

2.20 

8t 

I%J 

2.30 

2.27 

1.59 

3.00 

e 

(degl 

.648 

.578 

.483 

2.00 
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Hind Gust Rejection 


Response to s Tuo-Sigms Hoi l2onta1 Gust Impulse 
STOL - Method I and II 


outputs 
OLC ? 

3 ! 

1 

no 

5 

no 

3 

yes 

desinn 

criteria 

B 

max 

t ft/sec ] 

4.6 

■9QB 

mm 

6.75 

B 

max 

Iftl 

2.67 

2.43 

2.08 

20.0 

H 

max 

I f t/seo 1 

1.17 


.997 

6.50 

U« 

residual 

1 f t/seci 

-.001 

-.003 

-.032 

.675 

d 

residua] 

Iftl 

.015 

.019 

.075 

2.00 

• 

d 

residual 

I ft/sec 1 

-.010 

-.011 

-.008 

.650 


max 

IH] 

-.B84 

-.861 

-1.11 

6.00 

e 

max 

[deg ] 

-.623 

-.588 

-.443 

4.00 
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TABLE 4 

Mind Gust Rejection 


Response to a Tuo-Sigma Vertical Gust Impulse 
STOL - Method I and II 



outputs 


3 

s 

3 

design 


DLC ? 


no 

no 

yes 

criteria 


max 

1 ft/sec 1 

.577 

.536 

.476 

6.75 

d 

max 

Iftl 

2.51 

2.62 

1.68 

20. 

d 

max 

1 ft/sec 1 

1,03 

1.10 

.829 

6.50 

Ua 

residua] 

1 ft/sec 1 

.134 

.137 

.031 

.675 

d 

residual 

Iftl 

.528 

.519 

.090 

2.00 

• 

d 

residual 

I ft/sec 1 

-.055 

-.054 

-.017 

.650 

Si 

max 

ij:i 

-1.61 

-1.63 

-1.37 

6.00 

e 

max 

[deg ] 

-.455 

-.405 

-.356 

4.00 
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Figure 8: Response to Vertical Wind Gust (STOL - 3 sensors) 
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TABLE 5 

Response to Initial Offset in Airspeed (6.75 ft/sec) 
STOL - Method 1 and 11 


outputs 
DLC ? 


a max [ft/sec] 

max Ift) 

max [ft/sec] 


Ua residual Ift/sec] 

d residual [ft] 

d residual tft/ssc] 


St max 
6 max 


r/,1 

[deg ] 


3 

5 

3 

design 

no 

no 

yes 

criteria 

6.75 

6.75 

6.75 

6.75 

8.83 

8.29 

5.28 

20.0 

3.14 

2.83 

2.13 

6.50 

.555 

.584 

.133 

.675 

2.18 

2.21 

.391 

2.00 

‘.229 

-.229 

-.073 

.650 

-6.30 

-6.19 

-4.98 

6.00 

2.07 

1.65 

2.23 

4.00 
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TABLE 6 

Response to Initial Offset in Vertical Position (20 ft) 
STOL - Method I and II 


outputs 
OLD ? 

3 

no 

5 

no 

3 

yes 

design 

criteria 

Ua 

max 

If t/secl 

1.81 

1.64 

1,37 

6.75 

d 

max 

Iftl 

20 . 

20. 

20. 

20. 

• 

d 

max 

If t/secl 

-2.87 

-2.64 

-3.81 

6.50 

u. 

residual 

If t/secl 

.430 

.436 

.102 

.675 

d 

residual 

Iftl 

1.67 

1.65 

,298 

2.00 

• 

d 

residual 

1 ft/sec 1 

-.176 

-.172 

-.056 

.650 

St 

max 

[y,] 

-5. 10 

-5.16 

-2.57 

6.00 

e 

max 

1 deg 1 

-2.83 

-2.44 

-2.57 

4.00 
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6.3.2 Five Sensors Configuration 

Adding vertical acceleration (d) as the fourth sensor was not 
sufficient to meet the work-load criterion; and a fifth sensori axial 
acceleration (u) was added as well. Only the two additional output 
weighting and sensor variances were iterated to obtain the design 
examples presented below. The objective of these iterations uas to 
improve the work-load metric only, and no attempt uas made to improve 
the performance level. Varying the u weighting and sensor noise 
variance provided a much better control over the magnitude curve slope 
for the throttle channel than was available in the three sensor 

it 

configuration. The d weighting uas varied mostly to counteract the 
effect of the u weighting on the slope of the pitch attitude magnitude 
curve. 

Mind gust rejection and glide-slope capture performance are shown in 
the second column of Tables 2 to 6. As expected insufficient throttle 
effectiveness is still the limiting factor, and the improvement in 
performance is relatively small (5J{-10J{). As a result, the 
corresponding time history plots are not included, un the other hand. 
Figure 11 shows clearly the improvement in the work-load metric. The 
frequency response magnitude curve slope is -19.9 dD/decade (between 0.4 
and 4.0 rad/sec) for the pitch attitude channel. The corresponding slope 
for the throttle channel is -23.7 dB/decade, quite close to the desired 
value of -20. dB/decade. 



6.3.3 iIL£ Configuration (with 1 sensors) 
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The results presented so far, especially for the glide-slope capture, 
clearly indicate that performance is limited by the available throttle 
authority. Consequently, the use of the DLC mode where the spoilers are 
linked to the throttle was investigated. The spoilers to throttle ratio 
was chosen to effectively double the throttle authority (c.f.. Appendix 
A for the values of Z«p versus and Xgp versus X«}. The spoiler 
actuators are also very fast, thereby rartially compensating for the 
more sluggish thrust response of the engine. The loss of maximum usable 
lift is acceptable in view of the substantial stall margin of the CtSRA. 
It should be noted that the DLC mode is part of the original design of 
the aircraft, and its usefulness has been proven in flight. In 
particular, during carrier trials of the QSRA by the Navy ll7], the DLC 
mode was instrumental in bringing the pilot's uork-load down to 
acceptable levels (no flight director was used). 

The improvement in performance is quite dramatic, as shown in the 
third column of Tables 2 to 6. The corresponding time history plots are 
presented in Figurs 12 to 15. The work-load metric is shown In Figure 
11. The average slopes of the frequency response magnitude curves are 
-19.7 and -21.5 dB/decade for the pitch attitude and the throttle 
channels respectively (between 0.4 and 4.0 rad/sec). 

Last but not least, the above results were obtained using only three 
sensors, thereby simplifying the flight director implementation in the 
aircraft. This indicates that the faster acting spoilers, used in the 
DLC mode, reduce the need for lead equalization so that sensors to 
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measure higher derivatives of the aircraft's states are not required. 
As all design criteria were improved, no further iterations uere 
performed. The cost function weighting matrices are the same as those 
used in the three sensors configuration presented in Section 6.2.1. 
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6.4 C/STOL CONFIGURATION 

The C/STOL configuration (c,f.» Section 3.3) represents the 
transition between the conventional and the STOL configurations. The 
USB flaps are partially deflected to 30 degrees. Following the trade-off 
study results, DLC is used with three aircraft sensors. As this 
configuration is not as highly coupled as the STOL one* the value chosen 
for the DLC gain is only half of that used previously. Thus* a two 
percent increase in throttle will result in a one degree decrease in 
spoiler deflection. 


The results, presented in Tables 7 and 6. and in Figures 16 to 18, 
are similar to those of the STOL configuration with DLC. 


TABLE 7 

Wind Gust Rejection 


RMS Outputs and Controls 
C/STOL Configuration 


method 

I 

11 

Ua 

t ft/sec 1 

2.15 

2. 13 

d 

[ftl 

3,01 

3.81 

• 

d 

I ft/sec 1 

.940 

.872 

6t 

iJtl 

1.54 

1.53 

e 

[deg ] 

.563 

.559 
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TABLE 8 

Mind Gust Rejection and Glide-Slope Capture 
C/STOL - Method I and II 


criterion 

Uh 

Hh 

Uo 

Mo 

B 

max 

t ft/sec 1 

HI 

.688 

6.75 

1.88 

B 

max 

tftl 

BH 

1.68 

5.17 

20.0 

i 

mcx 

1 ft/sec 1 

m 

.859 

2.14 

-4.79 


residual 

1 ft/sec 1 

-.075 

.084 

.306 

.190 

d 

resi dual 

Iftl 

-.113 

.156 

.567 

.352 

• 

d 

residual 

[ft/sec 1 

.006 

-.024 

-.087 

-.054 

St 

max 

l%] 

-.879 

-1.38 

-4.73 

-3.82 

0 

max 

Idegl 

-.415 

-.616 

1 . 66 

-3. 52 


where 


Um ! response to horizontal gust impulse 

Um 5 response to vertical gust impulse 

Uo : response to initial offset in airspeed 

Mo s response to initial offset in vertical position 
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Chapter VII 

CONCLUSIONS AND RECOMtlENDATIONS FOR FUTURE RESEARCH 


7. 1 CONCLUSIONS 

Two methods for flight director control lau design have been 
presented. The first design method* proposed by Lev:nson* is based 
directly on the OCM concept. The second design method, proposed here, 
uses a fixed dynamic model of the pilot in a state space formulation 
similar to that of the OCM. The second design method also includes an 
explicit work-load metric. 

The several design examples presented in Chapter VI, clearly show the 
strong impact of throttle effectiveness on the performance and pilot 
work-load associated with manual control of a powered-lift STOL aircraft 
during approach. Use of DLC to increase throttle effectiveness greatly 
improves performance and reduces pilot work-load. 

Flight path reaulation in the presence of turbulence, initial capture 
of the glide-slope, and acceptable response to horizontal windshear, can 
be achieved, equally well, by both methods. 

The two design methods have the following features in commons 

1. They are based on a highly structured state space pilot/vehicle 
model, and use a Kalman-f il ter/LQ regulator. Their use requires 
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less expertise and engineering judgment than classical manual 
control design methods. 


2. The special requirements of the OCfl concerning pilot observation 
and motor noise-to-signal ratios can be easily achieved. 

3. Only three aircraft sensors were required to achieve the desired 
performance . 


Design method II has the follouing advantages over method 1 : 

1. An explicit work-load metric is included# which should improve 

the suitability of the flight director for manual control tasks. 
This metric is affected by the throttle effectiveness* and by the 
number of sensors used. With DLC only three sensors were 
necessary; without DLC* two additional sensors were required to 
ac>’eve the desired work-load level. Thus* a minimum sensor 
complement can be determined early in the design process which is 
not the case for method I. Iterating the weighing matrices and 
sensor noise levels to meet the work-load criterion has only a 
small effect on the flight director performance. Hence, the 

flight director can be designed in two consecutive phases : (1) 

performance improvement, and (2) work-load reduction. 

2. A simpler model of the pilot's observation and scanning process 
is required. The complexity of the observation model increases 
with the number of controls for the second design method whereas 
it increases with the number of outputs for the first method. 
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Also* the pilot's observation model used in method I is not 
relevant to the actual control task. The resulting control taus 
may have to be revised to account for the difference between the 
pilot observation noise levels used in the design phase and the 
actual sensor noise levels. 

3. The second method is more flexible in tnat any linear model of 
the pilot can be used, whereas the first method is limited to the 
OCn. Improved models of the human operator can be incorporated in 
the future with only a feu modifications. 


7.2 RECOmWENDATIOHS FOR FUTURE RESEARCH 

Design methods I and II should be evaluated in a piloted simulation. 

Of particular interest are : 

1. Pilot opinion ratings for the two methods. 

2. The performance levels achieved in the simulator as compared to 
those predicted by the two design methods. 

3. The correlation between the work-load metric proposed for method 
II and pilot opinion ratings. This could be done by comparing 
pilot opinion ratings for the three design examples presented in 
Chapter VI. or by generating additional examples with varying 
levels tf work-load. In addition to validating the proposed 
metric, xl-'ese data could be used as guidelines for simplifying 
the flight director control laws. 
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4. The robustness of the design to variations in the pilot model 
parameters : time delay# neuromuscular time constant* and pilot 

observation no'se levels. 
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Appendix A 

MODEL PARAMETERS FOR THE QSRA 


configuration j 

STOL 

C/STOL 

CTOL 

Uo 

Iktl 

70. 


130. 

h 

Iftl 

500. 


1000. 

Yo 

tdeg] 

-6. 

-6. 

0. 

6u*b 

Ideg ] 

50. 

30. 

0. 

Xu 

Isec* ’ 1 

-.074 

-. 066 

-.06 

Xm 

1 sec* ’ 1 

.131 

. 1056 

.094 

X. 

1 f t/sec^/Jil 

.031 

.367 

.22 

Xsp 

( f t/sec*/deg 1 

.034 

.012 


Zu 

Isec* ’ 1 

-.365 

-.304 

-.21 

Zm 

Isec* ’ 1 

-.485 

-.547 

-.862 

Z. 

I f t/sec^/JJJ 

-.327 

-.293 

-.22 

Zsp 

I f t/sec*/deg 1 

-.298 

-.506 


1/Tu 

I sec* ’ 1 

.3 

.237 

.342 

1/T„ 

[sec* ' ] 

.73 

.532 

.768 

Qu 

1 f t^sec* ’ 1 

35. 

45. 

31. 

Qn 

I f t*sec* ’ J 

14. 

19.9 

13.7 
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